Theoretical description of deformed proton emitters: 
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The newly developed nonadiabatic method based on the coupled-channel Schrodinger equation with 
Gamow states is used to study the phenomenon of proton radioactivity. The new method, adopting 
the weak coupling regime of the particle-plus-rotor model, allows for the inclusion of excitations in 
the daughter nucleus. This can lead to rather different predictions for lifetimes and branching ratios 
as compared to the standard adiabatic approximation corresponding to the strong coupling scheme. 
Calculations are performed for several experimentally seen, non-spherical nuclei beyond the proton 
dripline. By comparing theory and experiment, we are able to characterize the angular momentum 
content of the observed narrow resonance. 

PACS number(s): 23.50.-Fz, 24.10.Eq, 21.10.Tg, 21.10.Re, 27.60.-Fj 



OO 



> 

o 
l> 
o 
o 
o 



X 



I. INTRODUCTION 

Nuclei beyond the proton dripline are unstable against 
proton emission. Although formally unbound, some of 
these systems have rather long lifetimes, ranging from 
microseconds to seconds, due to the confining effect of 
the Coulomb barrier 

The past few years have seen an explosion of excit- 
ing discoveries in this field including new ground-state 
and isomeric proton emitters p-pland the first evidence 
for fine structure in proton decay . The focus of recent 
investigations has been on well-deformed nuclei which ex- 
hibit collective motion. These are of particular interest 
due to the interplay between proton emission and angular 
momentum. 

The theoretical description of long-lived proton emit- 
ters requires a detailed understanding of narrow reso- 
nances. Although proton radioactivity is a complicated 
A-hody phenomenon, much insight may be gained by 
considering the simplified problem of a single proton pen- 
etrating the Coulomb barrier of the core consisting of 
the remaining A-\ nucleons. It has been found that this 
simple one-body picture works surprisingly well. In many 
cases one has been able to determine the angular momen- 
tum content of the resonance and the associated spectro- 
scopic factor 

For spherical systems, there are many methods on the 
market which give similarly precise descriptions and, in 
many cases, one has been able to determine the angular 
momentum content of the resonance and the associated 
spectroscopic factor 

The array of theoretical tools available for deformed 
emitters is not as well developed. The existing ones 
fall into three general categories. The first family of 
calculations j^J^,^ is based on the reaction-theoretical 
framework of Kadmenskii and collaborators [|lO|. The 
second suite uses the theory of Gamow (resonance) 



states ||,|ll|-|T^. Finally, an approach, based on the time- 
dependent Schrodinger equation, has been introduced in 
Ref. 0. 

In all of these previous attempts, the strong coupling 
approximation of the particle-plus-rotor model has been 
used. The core is taken to be a perfect rotor with in- 
finite moment of inertia. This has the effect of (i) col- 
lapsing the rotational spectrum of the daughter nucleus 
to the ground state and (ii) neglecting the Coriolis cou- 
pling. Recently we have introduced a technique based on 
the weak coupling scheme which is free from these defi- 
ciencies ||l5[. Within this method, partial proton widths 
from different states of the parent nucleus to various fi- 
nal states in the daughter system can be calculated in a 
straightforward and consistent manner. 

We will begin in Sec. || by laying the theoretical frame- 
work for this work. Section III discusses the numerical 



methods adopted in our work. Section IV presents appli- 
cation of the method to the structure of deformed proton 
emitters. A critical analysis of the adiabatic and nonadia- 
batic methods is contained in Sec. Finally, conclusions 
are given in Sec. VI. 



II. THEORETICAL BASIS 

From a theoretical point of view, proton radioactivity 
is an excellent example of three-dimensional, quantum- 
mechanical tunneling. As such, the understanding of pro- 
ton emission is really a test of our knowledge of very nar- 
row resonances. Since the lifetimes which can be seen 
experimentally range from microseconds to seconds, the 
corresponding widths are extremely small; they vary be- 
tween 10"^^ MeV and 10"^^ MeV. Theoretical descrip- 
tion of such small widths requires high numerical accu- 
racy. In the following, the coupled-channel Schrodinger 



1 



equation method with Gamow states is outHned, and the 
proton-plus-core Hamiltonian is defined. 

A. Coupled-channel Equations 

The parent nucleus is described by the core-plus- 
proton Hamiltonian, 

H = Hd + Hp + V, (1) 

where Hd is the Hamiltonian of the daughter nucleus, Hp 
is that of the proton, and V is the proton-daughter inter- 
action. In the weak coupling scheme, the wave function 
of the parent nucleus is written as 

*,7M - (^'.7. ® *^J,/M ■ (2) 

This wave function is labeled by parity, total angular 
momentum J, and its projection M. In Eq. (||), u;^(r) 
[where a = (Jdlpjp) completely labels the channel quan- 
tum numbers] is the cluster radial wave function repre- 
senting the relative radial motion of the proton and the 
core, and yi j m is the orbital-spin wave function of the 
proton. The daughter wave function, ^j^m^, satisfies 

Hd<S>j,M,=Ej^<i>j^M^. (3) 

In the present formalism, the daughter's spectrum does 
not have to be known explicitly. Where possible, the 
energies Ej^ are taken from experiment; otherwise, the 
spectrum is modeled theoretically. Figure |l| shows a 
schematic diagram illustrating the energetics of proton 
emission from a J'^ state of an odd-Z parent nucleus to 
the ground-state rotational band of the deformed daugh- 
ter nucleus. 

As usual, the coupled-channel equations are obtained 
by inserting Eq. into the Schrodinger equation and 
integrating over all coordinates, save the radial variable 
r §|6|: 

+ E K.^"'M"f'M = o. (4) 

a' 

In Eq. (Q), Vaa is the diagonal part of the proton-core 
potential, Qj^ is the energy of the emitted proton leaving 
the daughter nucleus in the state Jd, and V^^, are the 
off-diagonal coupling terms. The Qj^ values follow from 
the spectrum of the daughter nucleus, Qj^ = Qo — Ej^, 
where Qo is the Qp value for the decay to the O"*" ground 
state (see Fig. 0). 

To illuminate the dynamics of the system, one can ex- 
pand the proton-daughter potential in multipoles [p^ , 

A 




FIG. 1. Schematic diagram illustrating the energetics of 
proton decay. A state of an oAA-Z parent nucleus (left) 
decays to those members of the ground-state rotational band 
of the even-even deformed daughter nucleus (right) which are 
in the Q window. The band members have excitation energies 
Ej^ relative to the ground state, and Qo is the Qp value for 
the decay to the ground state. As shown here, usually only 
a few channels are energetically open. The corresponding 
partial widths are indicated by Vj^ . 

The matrix elements ^, (r) can then be written in the 
simple, yet generic, form 

VlAr) = E^^W ^■M\Mx\\J'd) A{lpjpJd, I'pj'pJ'd, AJ). 

A 

(6) 

The factor A is purely geometric and comes from the 
proper coupling of angular momentum vectors. The re- 
duced matrix elements oi Mx contain all of the dynamics 
of the core. Since we consider only rotational nuclei in 
this paper, they are given by a simple expression fl^ 

{JdWMxWJ'd) - V^2 + 1 {J'dX KO\JdK). (7) 

To consider other excitation modes in the daughter sys- 
tem, one needs only change these reduced matrix ele- 
ments. 

To be a resonant state, the cluster radial wave func- 
tion must vanish at the origin and behave as an outgoing 
Coulomb wave, Oi = Gi + iFi, beyond the range of the 
nuclear interaction and the off-diagonal Coulomb inter- 
action, 

= Gip (VJa : rkj^ ) + iEi^ iVJ^ > rkj^ ) > (8) 

where fcj^^ = 2^Qj^/Tt? and rjj^kj^ — iiZe^/h^. These 
two conditions are only satisfied for a discrete set of 
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FIG. 2. Asymptotic behavior of three-channel wave func- 
tions corresponding to the J'^=3/2~ neutron resonance in 
"^Ho at energy £ = (0.378 - i0.0732) MeV calculated at 
/32=0.20. The solid line denotes the <8> 0"*" channel. No- 
tice the increasing envelope for increasing r. The dashed line 
labels the /5/2 <E> 2"*" channel function. The dotted line corre- 
sponds to the closed channel Z19/2 ® S"*". The wave function 
decays exponentially since Q < 0. The magnitude has been 
scaled so all three wave functions could be shown. 



complex wave numbers k. The generalized eigenval- 
ues of Eg. M) correspond to the poles of the scattering 
matrix |17|Jl8|] . The corresponding solutions are either 
bound or antibound states, £ = Eb < 0, with negative 
real energies and imaginary wave numbers k = ij > 
for bound and 7 < for antibound states ) , or resonance 
states, £ = Q — i^, with a nonzero imaginary part F ^ 0, 
and k — K — ij. 

The asymptotic behavior of these solutions is deter- 
mined by k; at a very large distance the outgoing so- 



lution is proportional to e . For resonance states, 
gifcr _ giK;rg7r^ j ^^iQ wavc function diverges exponen- 
tially. As discussed in Refs. [ pT|p^ , this seemingly un- 
physical feature of Gamow wave functions has a natural 
explanation in the fact that Gamow states do not rep- 
resent time-dependent wave packets but static sources. 
To illustrate the asymptotic behavior of Gamow wave 
functions. Fig. |^ shows three-channel wave functions cor- 
responding to a broad neutron resonance. 

Due to the divergent behavior at large r, one must de- 
fine a new normalization scheme for the Gamow states. 
Berggren proposed a new completeness relation, which 
includes Gamow states [0, by generalizing the scalar 
product. He introduced a bilinear basis set and a reg- 
ularization procedure {TZeg). With this generalization, 
the norm is 



J2neg [ui{r)]'dr = l. 



(9) 



A convenient method for regularization is to rotate r into 



the first quadrant of the complex r-plane beyond a cer- 
tain distance rmax- This is often referred to as the ex- 
terior complex scaling method. (For other regularization 
techniques, see Ref. Q.) 

Once we know the resonance energy and radial wave 
functions, there are several methods to calculate the 
width of the state. The most straightforward method 
is to take twice the negative of the imaginary part of the 
resonance energy. However, for the narrow resonances 
associated with proton emitters, the numerical accuracy 
needed to calculate Xm[£] is difficult to achieve. There- 
fore, approximate methods are often used. 

One possibility is to calculate the partial width for each 
channel from the current expression , 



<*(r)u„(r) - <(rX(r) 



where the sum of the partial widths, 

r(r-) = 5]r„(r). 



(10) 



(11) 



gives the total decay width (see Fig. ||) . Although values 
of Ta(r) depend on r in the region r<ras where the cou- 
pling potential terms are not negligible, T(r) is strictly 
independent of r by construction (T(r) — ~22m[£] = F) 
which reflects the flux conservation (continuity equation). 
Beyond the asymptotic radius, ras, the partial widths, 
Fa(r) have a negligible dependence on radius. We take 
^as ~ 40fm. Numerically, F(r) varies little with distance 
and differs by less than 0.1% from the F obtained from 
the imaginary part of the eigenvalue. 

The Gamow boundary condition given by Eq. (^) is 
usually written in the form 



"^ct (^as) 



k 



Jd 



OliVJd,ri,skj^) 
0/p(?7j^,rasA:.7j' 



(12) 



where ras is the channel radius. (The off-diagonal cou- 
plings are negligible beyond it.) Using Eq. (p^, the par- 
tial decay widths can be written at the point ras as 



EcK (^as ) — ^ 



\Uairas)\" 



|0/,(ry,;„A:,;,ras)PEaJo (r')Pdr' 

''as kj^ 



(13) 



If we neglect the very small imaginary part of kj^ , the 
square bracket in Eq. ( p^ ) is equal to —2i. Further- 
more, if we assume that for a very narrow resonance the 
imaginary part of Ua is very small (hence the general- 
ized normalization condition (^ is roughly equivalent to 
the "normal" normalization X]q' /o^"" \ua'[f')\^dr' « 1), 
then we end up with the approximate expression for the 
partial decay width: 
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Ta (^as ) 



\Ua{rss)['' 



(14) 



It is to be noted that Eq. ( [1^ ) and its approximate form 
( |r^ ) are valid only at the point r^s- The expression ([l^ ) 
was used in papers p^ , p^ , pOpl| ] . We emphasize that if 
the coupled equations are solved with the Gamow bound- 
ary condition, then the total width can be calculated at 
any intermediate point using Eqs. ( [l0| ) and The ex- 
pression ( p^ ) is very similar to that of the R-matrix the- 
ory (see below), but it relies on different approximations 
and boundary conditions than the R-matrix formalism. 

In the R-matrix theory, we also have a set of radial 
functions, ga- These functions are regular at the origin 
and satisfy the coupled equations but with the following 
boundary conditions: 



ffa(^as) _ „ 
'as / \ 



(15) 



where Ba are arbitrary real numbers. Due to the 
real boundary condition, the R-matrix eigenvalues are 
real numbers. In the R-matrix theory, the wave func- 
tion is normalized inside the sphere of radius ras, i.e., 
See /o"' \9a{'r)\'^dr = 1. Thomas has shown |22| that in a 
one- level approximation with appropriately chosen , in 
which the level shift is ignored, the position of the Gamow 
resonance corresponds to the R-matrix eigenvalue, and 
the width of the state can be calculated in the form given 
by Eq. ( p^ ) in which Ua{r) is replaced with ga{r). This 
R-matrix approximation works fairly well [ p3| , p4[ for very 
narrow Gamow resonances corresponding to the known 
proton emitters. For large values of the channel radius 
Tas, expression (|lj) is generally within 2% of the val- 
ues calculated explicitly from Eq. (^) or obtained via the 
current expression (^^. a detailed comparison of the 
R-matrix theory and the Gamow formalism for proton 
emitters will be given in Ref. |25|. 

The nonadiabatic approach allows for a straightfor- 
ward calculation of branching ratios. The partial width 
corresponding to the decay to a core state Ja is given by 



{ij} 



(16) 



where Fj^/j = F^ is given by Eq. (|l0|). Once the total 
width is known, the half-life for proton emission is 

?iln2 



Ti 



F 



(17) 



The use of the weak coupling scheme represented by 
Eq. (|^) has several advantages. First, excitations of the 
core are included in a straightforward manner. This en- 
ables us to study the proton decay from the rotational 
bands of the parent nucleus to the ground-state rota- 
tional band of the daughter nucleus. Furthermore, since 
the formalism is based on the laboratory-system descrip- 
tion [Hamiltonian (|l|) is rotationally invariant and the 
wave function vE'jm conserves angular momentum], the 
Coriolis coupling is automatically included. 



B. Strong Coupling Limit 

A great simplification to Eq. (Q) occurs if one considers 
all of the rotational states in the daughter's ground-state 
band to be degenerate (i.e., Qj^ = Qp for all J^). This is 
the limit of strong coupling where the moment of inertia 
of the daughter is taken to infinity. It is also the adiabatic 
approximation of Refs. [l^j2^ . 

In this limit, the coupled-channel equations, (^), re- 
duce to those for the intrinsic (i.e., deformed) Nilsson 
orbital 



Wo - 2^ yi^j^n 



(18) 



where 



ujKj.i, = V2 {-l)''+'Y.^JdQ3pK\JK)uj^i^^^. (19) 
Jd 

In Eq. ( p^ n^K—J is the angular momentum projec- 
tion on the symmetry axis. As seen from Eq. (^|), the 
strongly coupled intrinsic state contains contributions 
from all the cluster wave functions corresponding to dif- 
ferent core states. Since, as discussed by Tamura [ p^ , 
there is no dynamic coupling between the angular mo- 
mentum of the proton and that of the daughter nucleus 
(the daughter nucleus is perfectly inert during the proton 
emission), there exist infinitely many solutions obtained 
by combining jp and Jd- Since the core states are de- 
generate, all the solutions with J > il are degenerate as 
well. 



C. Model Parameters 

In this work, we assume that the average single-particle 
potential is approximated by the sum of a Woods-Saxon 
(WS) potential, a spin-orbit term, and a Coulomb po- 
tential. The axially deformed WS potential is defined 
according to Ref. We employ the Chepurnov pa- 

rameterization p8[ ; it is in good agreement with the pro- 
ton single-particle energy levels given in the systematic 
study |2^. Aberg et al. ||] discussed the effect of the 
optical model parameters on spherical proton emitters. 
They concluded that the uncertainties in the parameters 
affect the half- lives by, at most, a factor of 3. For spher- 
ical proton emitters, they concluded that the Becchetti- 
Greenlees WS potential |3^, commonly used in spheri- 
cal calculations for proton emitters, was better than the 
universal parameter set ]3l| ] (excellent for the description 
of structure properties of deformed rare-earth nuclei ]29| ] 
but having too large a radius to give a quantitative de- 
scription of the tunneling rate). Since, for the description 
of spherical proton emitters, the nodal behavior of radial 
wave functions plays a minor role |^ , the actual order of 
spherical shells does not really matter. 
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FIG. 3. Comparison of deformed single-proton levels for Z = 63, A'" = 68 predicted in three WS parameterizations. The left 
panel shows the Nilsson diagram calculated with the Becchetti-Greenlees set which yields poor ordering of the single-particle 
levels but good radial properties. The right panel is obtained with the "universal" set which yields good ordering of Nilsson 
levels but poor radial properties of wave functions. The center panel was obtained with the Chepurnov parameter set used in 
this work. This latter parameterization gives a very reasonable compromise between the radial and spectroscopic properties. 



However, in the case of deformed proton emitters the 
situation is different. While the radial properties of the 
optical model potential are still important, the proper 
ordering of spherical shells becomes crucial since it affects 
the fragmentation of orbital angular momentum caused 
by deformation. In this context, as illustrated in Fig. ||, 
the Becchetti-Greenlees parameter set performs rather 
poorly, while the Chepurnov parameterization offers a 
compromise between good radial properties and proper 
level ordering. 

Since within any mean-field theory the resonance en- 
ergy cannot be predicted with sufficient accuracy, follow- 
ing Refs. I^js), the depth of the WS potential is adjusted 
to give the experimental Qq value. The deformed part 
of the spin-orbit interaction is neglected; we do not ex- 
pect this to have a noticeable effect on the results [ p2[ . 
The off-diagonal coupling in appears thanks to the 
non-spherical parts of the WS and Coulomb potentials. 

Great care was taken to ensure that enough channels 
were considered in solving Eq. for proper convergence 
in the eigenvalues. As seen in the lower panel of Fig. ^, 
expanding the WS in spherical multipoles to order 8 is 
sufficient for convergence. However, to be on the safe 
side, a value of Amax = 12 in Eq. (|^) was used in all 
calculations. The number of partial waves that were 
needed in the decomposition of the proton radial wave 
function varies from system to system depending mainly 
on the angular momentum of the proton state. In gen- 
eral, all partial waves with /< 10-13 are needed. Since, in 
the nonadiabatic approach, the maximum proton angular 
momentum and the maximum daughter spin considered 



are closely related, the above condition corresponds to 
(>/d)max = 10 which was used for all calculations (see up- 
per panel of Fig. Since the high-spin channels are 
energetically forbidden, their exact placement is of mi- 
nor importance. Only the energy of the 2^^ level and, 
occasionally, the 4+ level have a profound effect on the 
resonance energy and other observable quantities. (For 
more discussion concerning this point, see Sec. 0) 

It is important in this work that we have a good rep- 
resentation of the ground-state rotational band in the 
daughter nucleus. In a few of the systems studied in this 
work, such spectra exists for < 10, which is enough for 
adequate convergence. However, for the most highly de- 
formed systems, the spectroscopic information does not 
exist. For these nuclei, we parameterize the ground-state 
rotational band as Ejj= KJd{Jd + 1), where k is ad- 
justed to the E2+ energy. In ^^^Eu, where fine structure 
has been seen, the E2+ energy is known. In other cases, 
systematic trends must be used. In actual calculations, 
we have used the NpNn scheme to estimate E2+ ■ 

III. NUMERICAL IMPLEMENTATION 

For realistic potentials, the radial Schrodinger equa- 
tion cannot be solved analytically but must be integrated 
numerically. For spherical potentials, one deals with 
a single radial equation instead of the full set (^). In 
Ref. |34j the code GAMOW was introduced, which uses 
the Fox-Goodwin method for solving the radial equa- 
tion. A more powerful method, the piecewise pertur- 
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TABLE I. List of Qp values, 2"*" state energies, and de- 
formation parameters for nuclei investigated. E2+ energies 
without a reference were estimated using the NpNn scheme 
of Ref. S. 
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FIG. 4. 



Dependence of the resonance eigenstate on various 

3 



expansion parameters. Calculations were done for the [411j ^ 
level in i^^Eu at a deformation of [32 — 0.37. The upper panel 
shows the resonance energy (solid line with squares and left 
scale) and the lifetime (dashed line with triangles and right 
scale) as a function of the number of included states in the 
ground-state band of the daughter nucleus. The lower panel 
shows the same except as a function of the number of spheri- 
cal multipoles used in expanding the deformed single-particle 
potential. 



bation, is used for the same purpose in Ref. ||35[]. The 
main features are similar in the two codes. The total 
r domain of coupled-channel equations (^ is separated 
into two parts. The first segment hes along the real axis, 
Ii — [0,rniax]. The other interval extends along the com- 
plex ray, I2 = [r max, fas], where Tas is complex and far 
enough away that at r^^ the asymptotic series of the out- 
going Coulomb wave, Oi{r]j^,raskj^), is a good approxi- 
mation. For a resonant state, the second integration re- 
gion must be complex for our regularization scheme given 
by Eq. (^. The rotation angle of r in I2 should satisfy 
the condition 



TT - arg(fcjj > arg(ras) > - arg(fcjj 



(20) 



so that the solution converges along the complex ray. 

For axially deformed V, the set of coupled-channel 
equations (jj) must be solved numerically. The piecewise 
perturbation method Pq] has been generalized for the 
coupled-channels case ]3^. A large value of r^s is used, 
which is far enough away that the off-diagonal terms of 
the coupUng matrix vanish and the asymptotic series for 



the Coulomb functions are accurate. At this point the 
coupled-channels equations decouple. For an initial kj^ 
value, one has to calculate the components of a "left" 
solution, M^, which vanish at the origin. These are inte- 
grated outwards to a matching radius, r^, in region Ii. 
The components of the "right" solutions, u^, are inte- 
grated inwards from Tag along the complex ray l2- At 
^■max, the integration path turns along the real axis to 
the matching radius, Tm. All components of the "left" 
and "right" solutions are linear combinations of linearly 
independent solutions of Eq. (|^) in the corresponding r 
regions. The two solutions and their derivatives with re- 
spect to r should match at the matching radius and form 
a set of functions which are continuous in r. This condi- 
tion gives a homogeneous set of linear equations for the 
unknown expansion coefficients of and u^. Non-trivial 
solutions exist only for the generalized kj^ eigenvalues 
where the determinant of the set of linear equations is 
equal to zero. For the initial value of kj^, the determi- 
nant is not zero; however, it is possible to find the zero 
of the determinant by iteration, e.g., using the Newton- 
Raphson method. For the known proton emitters, the 
width of the resonance is so small that extremely high 
numerical accuracy is needed to calculate the general- 
ized complex energy eigenvalue £. We have found that 
extended precision arithmetic must be employed to cal- 
culate the imaginary part of £ accurately. The width 
calculated directly in this manner matches well with the 
current expression (Eoh. 



IV. APPLICATIONS OF THE METHOD 

This section contains applications of the formalism to 
measured deformed proton emitters. For an easy orien- 
tation. Figs. H and g show the proton Nilsson diagrams 
characteristic of Z«55 and Z«67 nuclei, respectively. In 
our theoretical analysis, all Nilsson levels close to the 
Fermi level were investigated. The potential depth was 
always adjusted at each deformation so as to reproduce 
the experimental Qp value. Table || lists the Qp values, 
the energy of the 2+ states, and deformation parameters 
for all the nuclei investigated in this paper. 
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FIG. 5. Single-proton levels representative of odd-Z 
rare-earth nuclei with Z~55 plotted as functions of the 
quadrupole deformation /32- The hexadecapole deformation 
l34 was assumed to be proportional to f32 to give both the 
spherical and ground-state deformations. The Nilsson orbitals 
studied for ^^^^Cs, and ^^^La are marked by thick solid 
lines. 
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FIG. 6. Same as in Fig. | except for ^^^Eu and "^Ho. 



A. Description of Rotational Bands Built Upon 
Deformed Resonances 

As has been previously mentioned, a significant benefit 
of working in the nonadiabatic formahsm is the proper 
treatment of the ground-state rotational band in the 
daughter nucleus. This makes it possible to easily calcu- 
late the fine structure in the proton emission. The pres- 
ence of the rotational band in the daughter nucleus also 
gives rise to rotational bands built upon J — Q band- 
heads in the parent nucleus. In a previous work 
we discussed a rotational band in ^^^Eu built upon the 
J = I"*" level associated with the [411] | Nilsson orbital. 
The spacing of the levels in the parent nucleus follow 
nicely the expected J( J -I- 1) spacing with the same mo- 
ment of inertia parameter as assumed for the daughter 
nucleus. Small deviations from the J{J -t- 1) spacing re- 
sult from the Coriolis coupling. 
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TABLE II. A'-decomposition of the calculated band mem- 
bers of the [411] I band in ^^^Eu. The A" / | components 
arise from the Coriolis coupling. 



To verify that the calculated band structure indeed 
belongs to the same intrinsic Nilsson configuration, one 
can inspect the ^-decomposition of each rotational level. 
This is done by using Eq. ( p^ to project the nonadia- 
batic wave functions onto adiabatic states with good K. 
For the J=\^ orbital in ^^^Eu, the A'-decomposition is 
shown in Table ||. It is seen that the Ar=| dominates, 
although there appear small admixtures of other K com- 
ponents due to the Coriolis coupling. Note the presence 
of the K—^ which is forbidden in the strong coupling 
limit. 

A very different picture arises for the J—^ band built 
upon the [411]^ Nilsson orbital in ^'^^Ho. Its low-lying 
band members, through J=|, are shown in Fig. ^ In this 
case, we do not see the development of a strongly coupled 
band as in ^^^Eu, but rather two nearly degenerate de- 
coupled signature partners. This comes about due to the 
large decoupling parameter for this orbital. Since ^''^Ho 
is well deformed, we can consider the Coriolis interaction 
as a perturbation in the strong coupling approximation. 
For a, K — ^ band, first-order perturbation theory gives 



1 
2J 



J(J+l)-- + aa{~y+HJ+-) 



(21) 



where ad is the decoupling parameter. For a non-zero 
decoupling parameter, the J + \ odd levels are shifted 
against the J + \ even levels with a degeneracy setting 
in for \ad\ — 1. From studies of well-deformed and super- 
deformed bands in odd-Z rare-earth nuclei, bands built 
on the [411] i level are known to have a decoupling pa- 
rameter near —1 This nicely explains our pre- 
dictions and gives yet another verification that the weak 
coupling formalism properly incorporates the Coriolis in- 
teraction. It needs to be noted that the branchings shown 
in Fig. correspond to the proton emission only. In re- 
ality, the low-lying levels in these bands rapidly decay 
by gamma radiation (F^ ^ Fp); that is, the lifetimes of 
these states are much shorter. 

It is interesting to look in detail at the make-up of 
the cluster radial wave function and the partial widths. 
In the partial wave decomposition, the dominant compo- 
nents are those of the originating spherical state. For ex- 
ample, in ^^^La, the [422] | Nilsson orbital originates from 
a (77/2 spherical state. At a deformation of /32 = 0.33, 
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FIG. 7. Rotational bands in ^"^^Eu (bottom; built upon 
the [411]3/2 Nilsson level) and ^*^Ho (top; built upon the 
[411]l/2 Nilsson level). For ^^^Eu, the strongly coupled ro- 
tational band is predicted. For ^"^^Ho, the two decoupled, 
almost degenerate, sequences are calculated. Proton lifetimes 
and strongest branching ratios are indicated. 



the wave function still contains 60% of (77/2 distributed 
between the 2+ and 4+ daughter states. However, due 
to deformation, other partial waves with j > | also 
contribute: ^3/2 (9.6%), 4/2(9.0%), 59/2(14.2%), and 
Hi/2 (3.5%). Coriolis coupling introduces the si/2 par- 
tial wave with an amplitude of 4.2%. Although the radial 
wave function is a combination of components having dif- 
ferent angular momentum, the decay branches are easy 
to understand. The total width is governed by the high 
penetrability of low-^ partial waves. In fact, 97% of the 
width of this resonance is in the c?3/2 ^ 0+ channel. The 
remaining part comes from the ^3/2 (8) 2+ (0.6%) and the 
Si/2 (g) 2+ (2.3%) channels. 

The majority of decays investigated in this work have 
small branching ratios, less than ten percent. However, a 
few have quite large branching ratios to 2"*" states, includ- 
ing the possible decay out of the [532] | Nilsson orbital in 
^^^Eu which is predicted in this work to have the branch- 
ing ratio of 52%. The circumstances that lead to such 
large branching ratios are worthy of investigation. The 
[532] I orbital originates from an /111/2 spherical orbital. 
At a deformation of /32 = 0.32, the [532] | orbital con- 



sists mainly of /iii/2(75%), /7/2(18%) and only 1.9% of 
/5/2. There is an additional 0.8% of the iiT-for bidden P3/2 
component. The decay to the ground state can proceed 
only via the /5/2 component. Meanwhile, the decay to 
the 2+ state proceeds mainly through the P3/2 and /7/2 
waves; the former due to the lower angular momentum 
and the latter due to the larger make-up in the total wave 
function. The combination of a low-lying excited state, 
a lower angular momentum channel, and suppressed am- 
plitude of the /s /2 wave leads to the very high branching 
ratio this state would exhibit. 



B. Branching Ratios 

The main impetus behind this work has been the re- 
cently observed fine structure in the proton decay of 
"'^^^Eu [Q. The nonadiabatic formalism offers great ad- 
vantages over the strong-coupling approximation in cal- 
culating fine structure. The proper placement of the 
daughter states are explicitly included and the channels 
are now labeled with the proton's orbital and total an- 
gular momentum, Ij, and the angular momentum of the 
daughter nucleus, Jd- In one fell swoop, both the lifetime 
and partial widths are calculated. 

As was shown previously in Refs. [^,|l5| , for large defor- 
mations our calculations show little sensitivity to P2 and 
P4. This is because the spherical decomposition of the 
corresponding Nilsson orbitals varies little in this regime, 
provided that there are no crossings between the Nilsson 
orbitals of interest. The uncertainty due to nuclear defor- 
mation is usually smaller than that due to experimental 
uncertainty in the proton energy. In the less-deformed 
cases, there is a greater dependence of /32 and (34. In the 
1/2+ level in ^°^I, shown in Fig. |^, we see the effect of a 
level crossing near /?2 ~ 0.05. (This effect has been noted 
earlier in Ref [|l3||.) 



Table IH shows predicted half-lives, theoretical spec- 
troscopic factors, and branching ratios. The spectro- 
scopic factors have been estimated in the independent- 
quasi-particle picture. Note that the l/(il + |) factor 
present in the strong-coupling approximation is no longer 
needed. Our predictions for ^^^Eu and the ground and 



isomeric states in ^ Ho are unchanged from |15|. The 



ground state of ^^^Eu is consistent with the [411]|^^ as- 
signment. This is the same conclusion as in Refs . [Fp ^] 
but differs from the assignment of [413] | of Ref. [ [l3[ . 

In ^''^Ho, the assignments are straightforward: [523] | 
for the ground state and [411]| for the isomeric state. 
These match the assignments of Refs. |3l|l^. In ^°^I we 
find agreement with the [420] ^ with a deformation near 
/32=0.10. This agrees with suggestions of Refs. ||,|l^. In 
^^•^Cs we see a large admixture of K=l/2 in the J— 3/2 
wave function. Therefore, the asymptotic Nilsson label- 
ing is inappropriate, and only the total ang ular momen- 
tum is used to label the state in Table HI. The two or- 



bitals near the Fermi level in ^^■^Cs correspond to a pseu- 
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Orbital 



T"l/2 



b.r. 



0.20 



FIG. 8. Predicted lifetimes in for the [422]3/2 and 
[420]l/2 orbitals. The experimental lifetime is 110(5)^s Issl. 



dospin doublet; hence, strong CorioHs mixing of these 
levels. In the newly discovered proton emitter ^^^La, the 
experimental lifetime is 20(5) ms ^2). It appears that 
the [422] I assignment is best with a lifetime of 100 ms 
at a deformation of /?2=0.30, /34=0.11. 

There is currently a proposal [Q to search for proton 
emission from ^"^^Tb. Being in the region between ^^^Eu 
and ^"^^Ho, this nucleus is expected to be well deformed 
with /32«0.28. Using the Grodzins formula [^Q, we 
estimate the energy of the 2+ state in ^■^^Gd to be 165 
keV. Figure H shows the expected half-life as a function of 
Qp . It is expected that for lifetimes longer than the limit 
marked by the grey line, beta-decay will dominate Bn]. 



109 J 


[420] i 
[422] 1 


0.99 
0.99 


94.8 /is 
7.86 ms 
110(5)/is lH 


0% 
0% 


"^Cs 


[420] i 
J = 3/2+ 


0.52 
0.56 


0.66 /is 
34.7 /is 
16.7(7)/is § 


0% 
0% 


"^La 


[420] i 
[422] 1 
[541] 1 


0.32 
0.33 
0.61 


1.27 ms 
103 ms 
293 ms 
20(5) ms ^ 


0% 
3% 
4% 




[411] 1 
[413] 1 
[532] 1 


0.71 
0.52 
0.48 


34.0 ms 
184 ms 
3.90 s 
17.8(19) ms 


39% 
7% 
52% 
24(5)% |] 




[411] 1 
[523] 1 


0.70 
0.84 


14.6 fis 
19.1 ms 
3.9(5) ms |] 


0.8% 

6% 


141mjjQ 


[411] i 
[523] 1 


0.70 
0.84 


3.3 /is 
4.6 ms 
8(3)ms 1 


1% 
9% 



TABLE 111. Table showing the various orbitals for each 
system investigated in this work. Except for the weakly de- 
formed systems of ^"^I and ^^''Cs, the deformation dependence 
is much weaker than the uncertainty due to the experimental 
Qp value. The theoretical spectroscopic factor, half-life, and 
branching ratio to the 2+ states are shown. Experimental 
results (where available) are shown in bold type. 



C. Theoretical Uncertainties 

It should be emphasized that our method contains no 
adjustable parameters; there are a few parameters which 
are set by experiment. These include Qp and the place- 
ment of the lowest few levels in the ground-state band of 
the daughter nucleus. Since the higher levels are energet- 
ically forbidden, even if they are needed in the calculation 
to ensure proper convergence, the half-lives and branch- 
ing ratios are fairly insensitive to their placement. We 
shall now discuss the sensitivity of the calculated half- 
lives and branching ratios to various quantities used in 
the calculations. For concreteness, we will focus on the 
[411] I level in ^^^Eu. All other levels studied show simi- 
lar sensitivities. 

The largest effect on the lifetime comes from the Qp 
value. The Qp value for ^^^Eu is currently taken as 950(7) 
keV ]0] . The uncertainty of 7 keV leads to an uncertainty 
in the calculated lifetime of —7.5/ -I- 9.8 ms. This is a 
difference of roughly —22/ -I- 30%. Since a change in the 
Qp value also affects the energies of excited states, the 
change in branching ratio is much smaller. For the [411]| 



orbital, the effect is ±1.3 %. 

On the other hand, the placement of the 2+ level has 
a smaller effect on the lifetime but greatly influences the 
branching ratio. Based on Ref. [Q, the 2+ level in ^'^"Sm 
is placed at 121(7) keV. This 7 keV uncertainty changes 
the lifetime by ±4.0 ms (±12%). For the branching ratio, 
the corresponding error is ±6.7%. 

In the nuclei with significant branching ratios, little to 
nothing is known about the level structure in the daugh- 
ter system; hence, we had to assume a perfect rotor to 
assign energies to the states above the 2+. To check for 
the sensitivity to this assumption, we repeated some cal- 
culations assuming Ej^ = n' Jd{Jd + 1) ^ B[Jd{Jd + 1)]^- 
The anharmonicity factor, _B, has typical values around 
k7200 w 100 eV g|. This introduces a 1.0 ms shorten- 
ing of the lifetime and a reduction of the branching ratio 
of 1.2%. Both are much smaller than the influence of 
the Qp value or E2+ ■ So as long as the proper Qp value 
is used along with a good estimate of the first excited 
state, the remaining part of the spectrum needs only to 
be reasonably placed. 

Additional uncertainties can arise from the optical 
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H 10 



0.7 0.9 1.1 

Qp (MeV) 



1.3 



10-' 

10"- 

_0" 10"-' 
lO""* 



10" - 



10"- 



Non-Adiabatic 



Adiabatic 



Pl/2 P3/2 ^5/2 ^7/2 



FIG. 9. Predicted half-lives for "^Tb as a function of pro- 
ton Qp value. The assumed deformation is = 0.28 and the 
estimated value of E2+ is 165 keV. This figure is meant to 
identify the regions of Qp and T1/2 to look for this yet-unseen 
proton emitter. For longer lifetimes, serious competition from 
beta decay is expected. Above the grey line, the branching 
ratio for proton decay is predicted to be less than 10% ijlq]. 



model potentials. As discussed in Sec. [I C, we believe 



that the Chepurnov parameterization is the best current 
compromise. (It is noted here that better agreement 
between theory and experiment could, in principle, be 
achieved by fitting the optical model parameters to the 
properties, including proton decay data, of these drip-line 
nuclei.) As discussed in Ref. the lifetime of spherical 
proton emitters depends weakly on the nuclear structure 
details. Reasonable variations in radius and diffuseness 
parameters affect the lifetimes by less than a factor of 
about three. 



V. ASSESSMENT OF THE ADIABATIC 
APPROXIMATION 

As discussed in Sec all previous work on deformed 
proton emitters have made the adiabatic approximation 
(AD) |,|,0,|jll]-|l|,|^,|l|. The use of the nonadiabatic 
formalism for proton emission was first used by us in the 
recent Ref. [|l^ . The power of the nonadiabatic approach 
is apparent in several areas. First, due to the fact that 
the wave function is in the laboratory frame, the Cori- 
olis coupling is implicitly included. This allows for the 
inclusion of all the partial waves with j < J in the pro- 
ton's wave function. In particular, the Coriolis coupling 
can admix states with smaller Z-values, and consequently 
lower centrifugal barriers, into the proton wave function. 
Low-Z components, however small, can substantially af- 
fect the lifetime. Figure shows this effect for the [532] | 
deformed resonance in -'^^^Eu. Note that the P3/2 partial 
wave contributes only 0.3% to the total wave function. 



FIG. 10. Comparison of partial widths and wave func- 
tion amplitudes obtained in the nonadiabatic (left) and adi- 
abatic (right) calculations for the [532] | deformed resonance 
in ^''^Eu. Only the lowest few partial waves are shown. The 
upper panels show the spherical amplitudes \ca\'^ . Notice the 
presence of the A'-forbidden p-wave components in the nona- 
diabatic approach. The lower panel shows the normalized 
partial widths 7q = Fa/Ftot. In the non-adiabatic case we 
have summed over all possible daughter states. 



yet accounts for 15% of the decay width. 

In order to calculate the branching ratio in the adia- 
batic approximation, some ansatz must be used. Firstly, 
the partial width to the ground state is approximated by 
the width for the partial wave that matches the initial 
state. 



■pad 
^0+ 



1 



(22) 



For the excited states, a weighted sum over the possible 
partial wave components is used: 



ad _ 2 Jd -t- 1 



i 7. — 



|(J,0jf7|f^l7)|2 r^'i(QjJ, (23) 



where the partial widths have been calculated with a Qp 
value adjusted to that of the jj" state. This reduces to 
Eq. (22) for Jd = 0. This procedure will be referred to 
as the adiabatic corrected method (ADC). 

As can be seen in Table 0, in some instances the adi- 
abatic and nonadiabatic predictions are very close, like 
for ^^^"Ho and ^^^La. However, in other cases the dif- 
ferences are striking, like the factor of five difference in 
^^•^Cs. In those systems where the agreement is good, 
there is no admixture of lower-Z partial waves in the nona- 
diabatic formalism. 

It is also worth noting that lifetimes calculated in the 
full adiabatic method are usually shorter by a factor of 
up to four as compared to the ADC method. However, 
for the [411] I orbital of ^■^^Eu, there is a factor of 56 
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Nucleus Orbital 



^ (AD) 



(ADC) 



b.r. 



117 



La 



Eu 
i«Ho 
"i'"Ho 



[420] i 
J = 3/2+ 
[422] I 

[411] I 
[523] i 
[411] 1 



9.85 X 10"'' 
1.98 X 10"'' 
20.5 X 10"^ 
868. X 10"'^ 
5.95 X 10-3 
1.94 X 10-*^ 



24.1 X 10"'' 
7.02 X lO"'' 
90.3 X 10^3 
48.7 X 10^3 
6.66 X 10^3 
3.43 X 10"^ 



0% 
0% 
0.3% 

37% 
3% 



TABLE IV. A calculatio n of t he proton lifetimes using the 
adiabatic formalism of Sec. 



IT E . All lifetimes are in seconds. 



The label (AD) corresponds to calculations in the adiabatic 
approximation, i.e., it includes all degenerate final states. The 
column ADC corresponds to the corrected adiabatic approxi- 
mation, in which it is assumed that the decay goes only to the 
ground state of the daughter nucleus. Eg. (|22|). The branching 
ratio is calculated using Eqs. ( p^ ) and (^3|). For the nonadi- 
abatic model predictions, cf. Table 



difference. This results from the large component 
in the corresponding Nilsson model function. 

In about half of the cases studied, the adiabatic ap- 
proximation, particularly with angular momentum con- 
servation enforced by hand, gives results similar to the 
nonadiabatic method. In the rest, the difference can be 
large. 



VI. CONCLUSIONS 

The state-of-the-art coupled-channel formalism has 
been extended to include excitation modes in the daugh- 
ter system. The weak-coupling scheme applied allows 
us to work in the laboratory reference frame. The exact 
treatment of excitation spectrum in the daughter nucleus 
also allows a consistent calculation of branching ratios. 

As could be expected, significant branching ratios are 
expected only for well-deformed nuclei where the first 
excited state of the daughter nucleus lies low in energy. 
The Coriolis mixing of states with lower orbital angular 
momentum can enhance the decay to the excited state, 
e.g., the decay of the [532] | orbital in ^'^^Eu where the 
branching ratio to the 2+ state is predicted to be as large 
as 52%. 

In the case of spherical proton emitters, the proton 
separation energy and orbital angular momentum have 
the largest effect on lifetimes The detailed nuclear 
structure plays a minor role. In the deformed case, the 
placement of the single-particle levels also has a signif- 
icant effect. This is true in both the adiabatic formal- 
ism, where mixing occurs with higher-^ states, and in the 
nonadiabatic formalism where the Coriolis coupling can 
mix states with lower-Z also. 

We have been able to calculate the placement and de- 
cay properties of excited levels in the ground-state band 
of the parent nucleus. As shown in Fig. 0, both strongly 
coupled and nearly degenerate decoupled bands are pre- 



dicted, depending on the nature of the band head. While 
we have calculated the proton decay half-lives of the 
excited band members, they are all much too slow to 
compete with in-band gamma-decay. In all cases investi- 
gated, by comparing theoretical predictions with experi- 
mental half-lives and branching ratios (where available), 
we have been able to identify the Nilsson orbital which 
the proton occupies. We have also been able to discern 
the angular momentum components of the proton wave 
function. 

While our calculations for well-deformed nuclei give 
a quantitative agreement with experiment, it would not 
be proper to apply the present model to vibrational or 
transitional nuclei such as ^^^La. For this, one needs to 
redefine the coupling matrix elements in Eq. (H). Calcu- 
lations along these lines are in progress. 
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